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Abstract: Krzakala, Ricci-Tersenghi and Zdeborova have shown recently that the ran- 
dom field Ising model with non-negative interactions and arbitrary external magnetic 
field on an arbitrary lattice does not have a static spin glass phase. In this paper we 
generalize the proof to a soft scalar spin version of the Ising model: the Ginzburg- 
Landau model with random magnetic field and random temperature-parameter. We do 
so by proving that the spin glass susceptibility cannot diverge unless the ferromagnetic 
susceptibility does. 

1. Introduction 

A widely studied class of disordered systems in statistical physics consist in adding 
random impurities to a field coupled with the order parameter. A textbook example of 
such a system is the random field Ising model (RFIM), introduced by [ 1 ], that has been 
a very useful playground for theoretical ideas. The Hamiltonian of the standard RFIM 
reads 



where all the non-zero interactions are ferromagnetic, i.e. Jy > 0. The N Ising spins 
Si — ±1, i — 1, . . . , N, are placed at the vertices of a graph (lattice), and hi are 
quenched random magnetic fields. The fact that all the interactions Jy are non-negative 
is fundamental, it means that in the absence of the fields there is no explicit frustration 
in the problem. 

The case where the graph of interactions is a finite-dimensional lattice and where 
the fields are taken from a Gaussian distribution with zero mean and a variance Hr 
has received a lot of attention. Of particular interest is the phase diagram in the T- 
Hr plane, where T is the temperature. Several authors have suggested, based on non- 
rigorous field theoretic arguments, that there exists an equilibrium spin glass phase in 
the three-dimensional RFIM, that is a phase with a random frozen ordering [2,3,4,5, 




(1.1) 
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|6|. These suggestions were disproved rigorously in [7 1 for the RFIM defined by Hamil- 
tonian dl.U . In particular [7| showed that for the RFIM d 1 - 1 b a special case of the 
Fortuin-Kasteleyn-Ginibre (FKG) inequality 1 8 1 implies that the spin glass susceptibil- 
ity is upper-bounded by the ferromagnetic susceptibility. Since the spin glass suscepti- 
bility diverges in the whole spin glass phase, a spin glass phase can not exist away from 
the ferromagnetic critical point/line in the RFIM. 

The field theoretic approach of [2 3 4 5 6], however, was not formulated with the 
Ising spin Hamiltonian ( II. U but instead with the soft-spin description of the random 
field model. This is the well-known Ginzburg-Landau model (or the so-called <j> 4 - 
theory) which is defined by the following Hamiltonian 

Hn = - ^2 Jij^i ~ z2 + /2 Ti< & + /2 Ui ^i ' ( L2 ) 

ij i i i 

where (pi are now real numbers, </>j G M, and the interactions are ferromagnetic, Jy > 
(this will be the case in the whole article). 

The generalized model i ll .21 includes several special cases. The Ising model is re- 
covered in the limit where r, = — 2ui and ui — > oo. The most common random field 
model is obtained when hi are random variables while r.i = r and m = u are fixed, 
and Jij = or = 1 depending if the spins ij interact or not. Another version that 
was considered in the literature, the random temperature model, is when are random 
while hi = 0, Ui — u and E {0, 1}. The existence of a spin-glass phase was also 
predicted in the random temperature model 191 1101 . based again on some non-rigorous 
arguments using perturbation theory; this result was, however, questioned in ifTTI . 

Our results work even for the slightly more general Hamiltonian 

Un = - £ J <i&& + £ MM ' (L3) 

ij i 

where Jy > and the local constraining potentials are arbitrary analytic functions, 
but for the requirement that the partition function 

/■oo n 

z N = n^ e_ ^ Wi}) ' ^ 

must exists for any N G N. This is the most minimalist requirement, since the non- 
convergence of the integral in d 1 -4b would make the Gibbs-Boltzmann measure ill de- 
fined and the model would not be a physical one. 

The Gibbs-Boltzmann average at temperature T = fi^ 1 is defined by 

{A) (N) = ^ TTd^e-^«^», (1.5) 

The superscript (TV) on the angular brackets will be written explicitly only when the 
size dependence is crucial, while the temperature dependence is always made implicit. 
Connected correlation functions are defined as 



(AB) C = (AB)-(A)(B) 



(1.6) 
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It is worth noticing that the convergence of the integral in ( 11.4b ensures that the 
partition function Zn is an analytic function of the coupling constants J,j for any finite 
value of N. Then the derivative 

exists as well for any pair of indices i, j, and this implies that single-variable marginal 
probability distributions have the first and the second moment, ((pi) and (<f>f). Actu- 
ally in soft-spin models used in the literature, such as the spherical model and the (p 4 
model, single-variable marginal probabilities decay exponentially fast for large values 
of (pi, and so all the moments (<p!?) exist. However, our proof only requires the first two 
moments to exist. 

Note also that any type of lattice can be encoded in model ( 11.3b by setting Jy = 
if spins i and j do not interact. 

The main contribution of this paper is a rigorous proof that the soft-spin random- 
field random-temperature model defined by ( 11.3b does not have a spin glass phase as 
long as the interactions are ferromagnetic (non-negative). This generalizes the result of 
0. 



2. Definitions of susceptibilities 

We define the ferromagnetic and the spin glass phases using the properties of the ferro- 
magnetic and spin glass susceptibilities. 

The order parameter that characterizes a ferromagnetic transition is the magnetiza- 
tion Til = X)i(^i)/-^- However, a non-zero magnetization does not imply a ferromag- 
netic phase. Indeed, m > even at large temperatures when a uniform positive external 
magnetic field is applied. A convenient way to characterize the ferromagnetic phase is 
to define the ferromagnetic susceptibility as 



N 

where 



(2.1) 



(2.2) 



VW) - W 2 

are the fluctuations with respect to the average values, normalized by the variances. 

In the thermodynamic limit (N — > oo), x% (°°) i s finite m the high temperature 
(T > T c ) paramagnetic phase and it diverges approaching the ferromagnetic critical 
point from above (T \ T c ). Right at the critical point (T = T c ), x°f( n ) diverges 
with N — > oo signaling that the system is critical, i.e. has long range correlations 
between fluctuations of its variables. Unfortunately, the ferromagnetic susceptibility 
Xf(N) diverges with N also in the whole low temperature (T < T c ) ferromagnetic 
phase: however this divergence is not due to criticality (i.e. long range correlation of 
fluctuations), but only because below T c two ferromagnetic states coexisfl 



1 In the presence of two or more equivalent states, an appropriately chosen perturbation, although of in- 
finitesimal strength, may induce a macroscopic change of state, thus leading to an infinite susceptibility. 
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Given that we are interested in finding critical points and critical lines where a phase 
transition takes place, we would like to measure an observable that diverges only at 
criticality, and so we consider the following ferromagnetic susceptibility 

Yp = lim lim Ypih, N) = lim lim — S^iSfa 8<hj) , (2.3) 

y 

where ft is an auxiliary uniform magnetic field (in practice one needs to add a term 
— ft J2i & m me Hamiltonian). Due to the order of the limits in (12.31 1. below T c , the 
infinitesimal external field ft makes the two ferromagnetic states no longer equivalent, 
and consequently \f is finite everywhere, but at the critical point T c (which is indeed 
defined as the point where xf diverges). 

In general to define a susceptibility that diverges only when a critical state is present 
one should explicitly break (by adding infinitesimal perturbations) all the symmetries of 
the Hamiltonian. In our case, the Hamiltonian (11.31 ) is very general, but the first term is 
invariant under the transformation fa — > — fa Vi. In case the potentials too are invariant 
under such a transformation, fi(4>) = /j (—</>), then the infinitesimal auxiliary uniform 
field in (12.3b is strictly required. 

The spin glass phase is characterized by a freezing of spins in random directions 
ifTZl . hence the spin glass susceptibility is defined as 

h\,0 N— >oo h\0N— »oo iV ^— ' 

ij 

Again we use the infinitesimal auxiliary external field to break the <j) — > —<fi symmetry, 
if present. The susceptibility xsg is related closely to what is measured in simulations 
and experiments lfl3l . and it is predicted to diverge at the critical point in spin glass 
theories such as replica symmetry breaking 04), or the droplet description [15|. In a 
spin glass phase \SG is infinite, because of the presence of at least two stated related by 
symmetries, which are not broken by the auxiliary field. For this reason we can define 
that a system is in a spin glass phase if and only if the ferromagnetic susceptibility (12.31 l 
is finite, while the spin glass susceptibility ( 12.4l i is infinite. 

More precisely the computation of these two susceptibilities must proceed by first 
taking the thermodynamic limit in the presence of the external field, XF(h, oo) and 
Xsg (ft) oo), and then studying the limit ft, \ of these two functions. If such a limit 
exists, then we say that the susceptibility is finite and we are away from the critical 
point, while if a divergence is found while decreasing ft, then we say that the suscepti- 
bility is infinite. 

In the next Section we prove that xsg (ft, N) < XfQi, N), for any value of ft and 
N, thus excluding the possibility of a spin glass phase (defined by xf < oo and xsg ~^ 
oo) in the model (11.31 1 in the absence of explicit frustration in the couplings. 

3. Results 

We start by proving a generalization of the 2nd Griffith's inequality 1 16 1. The follow- 
ing Lemma is also a consequence of much more general FKG inequalities ll MITl . we, 
however, find useful to present an independent and more elementary proof. 



2 The number of states depends on the model and for some models, like the 3D Edwards-Anderson model, 
it is still a matter of debate. 
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Lemma 3.1. In the model defined by the Hamiltonian ( I7.3D with non-negative cou- 
plings, Jij > under the Gibbs-Boltzmann measure e~P HN /Zn the correlation 
between fluctuations of any two variables is non-negative and bounded by 1, 

<(<%<%)<! Vi,j. (3.1) 



Proof (Lemma UTPb . Let us prove first the second inequality in ( 13.11 ). From the definition 
d2.2l > of the relative fluctuations we have that (Scf>f ) = 1 for any i. Moreover for any 
pair of indices i, j we have that 



< ((<% - 5^f) = (6$) + (Stf) - 2{8cj> i Stf>j) = 2(1 



bj)) (3.2) 



from which (5(f>i 8<fij) < 1 follows. 

In order to prove the first inequality in ( 13.11 ) we notice that it is equivalent to the 
inequality 

(Mi)c > , (3.3) 

thanks to the fact that all denominators in the definition ( 12. 2t of 5<j>i are positive and 
can be canceled without changing the sign of the correlation. 

Then we prove ( 13.31 ). by induction in the system size. In a system of TV = 1 spin 



A?>c W >0, 



(3.4) 



since the variance is always non-negative. Then we assume the property to hold in a 
system of N spins and we consider a system with N + 1 spins. The Hamiltonian of that 
system is related to the TV-spin system as 



N 



Wn+X — Wn — ^ JN+l,i4>N+l4>i + /jV+l(^JV+l) 



(3.5) 



We denote 



P(x) 



->N+1 J-oo 



JV 



(3.6) 



exp 



N 



-pH N {{4>i\) + JN+l,iX<Pi ~ PIn+i{x) 



Let us denote the thermodynamic average in a modified external magnetic field as 



(N) 



JUL 



(3.7) 



The connected correlation between spins 4>n+i and </>,; in the N + 1-spin system can 
then be rewritten as 



N+1 fr)i N+ V= 6xxP(x)(^- dyyP(y) dxP(x){^ 



— OO 

oc 



dx 



dyyP(y) 



P(x)(&>W 



(3.8) 
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We can then use the following inequality: For any real non-decreasing function g(x) 
such that 

/■OO 

(3.9) 



(3.10) 



Proof of this statement is elementary, function g(x) has to be non-positive for some 
x < xq and non-negative for x > xq. Since f(x) is non-decreasing one has 



da; g(x) = , 
and any non-decreasing function f(x) one has 

dxg{x)f(x) >0. 



da; 1^(^)1 /(^) < / dx \g(x)\ f(x) 



(3.11) 



from which ( 13.10b follows. We observe that (4>i)x is a non-decreasing function of x. 
Indeed 



d<&> 



(N) N 



dx 



(3.12) 



where the last inequality follows from the inductive assumption and since Jn+ij > 0. 
And that 



dx 



dyyP{y) 



P(x) = 0. 



Hence ( 13.101 ) implies 



(3.13) 



(3.14) 



We proceed similarly for the correlation function between two spins that were al- 
ready present in the TV-spin system 



dxP(*)(^)W 
dxP(x)(^)i N J 

1 

dxP(x)(<lH)W 



dxP(x)fc)W / dyP(y)(^ 



/ dyP(„)<^)W 



(3.15) 



where the first term is non-negative by the inductive assumption, and the second term 

is non-negative according to ( 13.101 ), because ((f>i) x and (4>j)x are non-decreasing 
functions of x and 



da; P(x) 



0. 



Hence 



(3.16) 
(3.17) 



Based on the previous Lemma we can now easily state the main result of the present 
paper. 



(Mi N+1) > o. 

This concludes the proof of Lemma [3TT1 □ 
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Theorem 3.1. In the model defined by the Hamiltonian ( I7.3D with non-negative cou- 
plings, Jij > Wi,j, under the Gibbs-Boltzmann measure e~^ N /Zff, the spin glass 
susceptibility Xsc{h, N) is always upper-bounded by the ferromagnetic susceptibility 
XF(h, N). Consequently the model does not posses a thermodynamic spin glass phase. 

Proof (Theorem \3.1h The hypothesis of the present Theorem are the same as those of 
Lemma [3~Tl with the only difference that to properly define the susceptibilities we need 
to add the external auxiliary field term to the original model Hamiltonian. Then Lemma 
I3.1l can be used only if 

r oo N 

Z N (h)= JJd^e-^^W+^E^ (3.18) 

exists also for h > 0, and this is easy to prove. Indeed Zn(0) exists (otherwise the 
Gibbs-Boltzmann measure would be ill-defined) and also the first two derivatives of 
Z^(h) with respect to h exist (because (<f>i) and (<j)f) exist): so Zjsr(h) can be continued 
in a neighborhood on h = 0, that we call Sq, and this is enough to take the limit 
h \ that is required to define properly the susceptibility. Please note that the region 
So coincide with K for all the models used in the literature, such as the spherical model 
and the cj) 4 model. 

Given that the hypothesis of Lemma [3~T1 are satisfied in Sq, we can make use of 
inequalities in ( 13.11 ) and find that 

(S^) 2 < (SfcSh) =^ XSG(h,N)< X F(h,N) (3.19) 

for any value of N and h 6 Sq. Even in the thermodynamic limit the inequality holds 

XSG(K oo) < XfQi, oo) (3.20) 

and so the spin glass susceptibility can not diverge if the ferromagnetic one stays finite. 

In other words, from the definitions given in the previous Section it is clear that 
if a thermodynamic spin glass phase exists, then for a sufficiently large value of N 
and a sufficiently small value of h the spin glass susceptibility must be larger than the 
ferromagnetic one and this would violate the inequality in (13.191 ). Then we conclude that 
a thermodynamic spin glass phase does not exists in the model defined in the hypothesis. 
□ 



4. Discussion 

We have shown rigorously that there is no spin glass phase in the scalar soft-spin 
random-field random-temperature Ginzburg-Landau model with ferromagnetic inter- 
actions defined by (11.31 ). This shows that with two-body interactions and a scalar order 
parameter one cannot obtain a genuine spin glass phase at equilibrium without explicit 
frustration in the couplings (another possibility to frustrate the system is to impose a 
non-equilibrium value of magnetization, see Q). 

Our proof contradicts the conclusions of some works that used field theoretic ar- 
guments [2,3,4,5,6,9,10]. It is yet to be discovered where the problem lies in those 
approaches. One possibility to consider is that the spin glass instability could be an 
artifact of some truncation in the perturbative expansion. For some of these works 
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the discrepancy may stem from the use of vectorial soft-spin models instead of scalar 
ones. Another possibility, that is related to what was suggested recently in ifTHl . is that 
the observed "replica symmetry breaking" instabilities arise only in disorder averaged 
quantities and never in the thermodynamic limit of a single instance quantities. These 
instabilities would then not be equivalent to the divergence of the spin glass suscep- 
tibility (which we prove impossible out of the ferromagnetic critical point), but they 
could instead be connected to some subtle non-self-averaging effects between different 
realizations of the system. Indeed all the above-mentioned works considered a "repli- 
cated" field theory, that is, a field theory averaged over many realizations of the disor- 
der. The divergences that they observed could hence be coming from strong sample to 
sample fluctuation. The fact that some non-self-averaging is present in the RFIM has 
been suggested by Parisi and Sourlas [19]. They argue that the correlation function, or 
equivalently the ferromagnetic susceptibility, of the RFIM is non self-averaging in the 
critical region, and they argue that this was the source of the problems with perturba- 
tive expansions. Note, however, that such simple non-self-averaging effects can only 
take place at the ferromagnetic critical point in any finite dimensional system. This is 
a straightforward consequence of a theorem by Wehr and Aizenman ||20| who proved 
that any extensive quantities (such as the ferromagnetic susceptibility away from the 
critical point) is self-averaging in finite dimensional systems. In other words, if this ef- 
fect was the one observed in the field theoretic approaches, it has to be limited to the 
ferromagnetic critical point itself. 

Finally, it would be very interesting to see if our proof can be generalized further. 
There are two interesting counter-examples that seem to put strong limits to such gen- 
eralizations. Matsuda and Nishimori (private communication) showed that a random 
field Ising model with 3-spins interactions on the Bethe lattice can have a spin glass 
phase. And so moving beyond pairwise interacting models seems impossible in full 
generality. Moreover Parisi (private communication) provided an interesting example 
of a pairwise interacting n = 2 component vector spin system where the two point 
connected correlation can be negative even if all couplings are positive. It is a chain of 
spins with an external field that smoothly rotates by 180 degrees along the chain, such 
that the field on the last spin is opposite to field on the first spin. If the field strength is 
strong enough, each spin will be mostly aligned along the local field and will thermally 
fluctuate around this position. However, given that the extremal spins are in opposite 
directions, their thermal fluctuations will be negatively correlated. This is a very spe- 
cific configuration which may not happen in typical samples, but its existence implies 
that the proof strategy presented in this paper cannot be straightforwardly generalized 
to vector spin models. 

Acknowledgment: We thank G. Parisi, H. Nishimori, F. Toninelli and for very useful 
comments and enlightening discussions. 
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